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Abstract 

The Telescope Array (TA) is an experiment to observe 
Ultra-High Energy Cosmic Rays (UHECRs). TA’s re¬ 
cent results, the energy spectrum and anisotropy based 
on the 6-year surface array data, and the primary com¬ 
position obtained from the shower maximum (Xmax) are 
reported. The spectrum demonstrates a clear dip and 
cutoff. The shape of the spectrum is well described 
by the energy loss of extra-galactic protons interacting 
with the cosmic microwave background (CMB). Above 
the cutoff, a medium-scale (20° radius) flux enhancement 
was observed near the Ursa-Major. A chance probabil¬ 
ity of creating this hotspot from the isotropic flux is 4.0 
cr. The measured <Xmax> is consistent with the pri¬ 
mary being proton or light nuclei for energies 10^^-^ eV 
- 10^^-^ eV. This report is to appear in the proceedings of 
the ISVHECRI-2014 symposium at CERN, August 18th- 
22nd, 2014. 

1 Introduction 

A cutoff of the cosmic ray energy spectrum around 
10^^ eV was suggested by Greisen, Zatsepin and Kuzmin 
(GZK) in 1966 HI just after the discovery of the cos¬ 
mic microwave background (CMB). When the energy of 
cosmic ray proton approaches to ultra-high energies, pion 
photo-production with the CMB starts. The creation and 
decay of the nucleon resonance deprives a significant en¬ 
ergy from the propagating proton. This process produces 
a cutoff and a horizon for UHECRs arriving at the Earth. 
Precise measurements of the pion photo-production cross 
section by accelerator experiments, and of the CMB radi¬ 
ation temperature by space experiments allowed a precise 
prediction for the cutoff energy and the horizon, assum¬ 
ing certain models for cosmic ray production and its cos¬ 
mological development Thus, a search for the GZK 
cutoff became an important subject for cosmic ray exper¬ 
iments, and a search for UHECR sources within the GZK 
horizon became a realistic hope. 

Searches for the GZK effect by ground based air shower 
arrays have been associated with the experimental diffi¬ 
culties such as the very low arrival rate, one event per year 
on an area of 100 km^ at ~10^^ eV, and large ambiguities 
in the determination of primary energy and composition 
by the sparse sampling of shower particles arriving on the 
ground. 

A new air fluorescence technique for measuring UHE¬ 
CRs was established by the Ely’s Eye (EE) and High Res¬ 
olution EE (HiRes) EIH by 2000. In this technique, the 


longitudinal development of air shower in the atmosphere 
was recorded by fast imaging telescopes, and the primary 
cosmic ray energy was determined by integrating the pro¬ 
file of the shower development. A stochastic measure¬ 
ment of the primary composition became possible by ob¬ 
serving the average and fluctuation of Xmax over many 
shower samples. 

The Telescope Array (TA) employs an array of Surface 
Detectors (SDs) covering a large ground area of 700 km^, 
together with a set of Eluorescence Detectors (EDs) over¬ 
looking the SD area ||5] |6l. It is installed in the West 
Desert of Utah, USA, and has been taking data since 2008. 
Whenever possible the SDs and EDs make simultaneous 
measurement of the event, however the operation of ED 
is limited to moonless clear nights and its duty factor is 
at a level of 10%. The SD runs with close to 100% duty, 
collects higher statistics, and its sky sampling is uniform 
in equatorial longitude. A similar type of hybrid experi¬ 
ment was installed at the Pierre Auger Observatory (PAO) 
in the Argentinian pampas 0, and has been taking data 
since 2004 with 3000 km^ SD array. 

The search for the GZK effect, and understanding the 
nature of UHECRs is rapidly proceeding with the two 
large-scale hybrid experiments, TA and PAO. We report 
here the recent results from TA as of summer 2014. 


2 TA Experiment 

The TA is located in Utah, USA at a latitude of 39.3° 
North and a longitude of 112.9° West. The altitude is ap¬ 
proximately 1400 m above the sea level. The TA encloses 
507 SDs, 3 stations of EDs and several calibration devices 
on site (see EigHJ. 

The TA SDs are deployed in a 1.2 km grid covering the 
ground area of 700 km^. Each SD is made of 2 layers 
of 1.2 cm thick, 3 m^ large plastic scintillaters overlaid 
on top of each other 0. Each layer is read out indepen¬ 
dently with a PMT using wavelength shifter fibers. The 
waveform is locally recorded by 12bit, 50MHz EADC for 
any hit with more than 0.3 Minimum Ionizing Particles 
(MIPs) in both scintillators in coincidence. These events 
are also locally histogrammed and used for the calibra¬ 
tion. The stored waveforms are read out using a 2.4 GHz 
wireless network when more than 3 MIPs are recorded in 

3 adjacent SDs within 8 /iS. 

The TA ED is composed of 38 fluorescence telescopes 
installed in 3 ED stations at Black Rock Mesa (BRM), 
Long Ridge (LR) and Middle Drum (MD) surrounding 
the SD array. The ED at MD station, ED(MD), is in¬ 
strumented with the refurbished HiRes telescopes with 
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Figure 1: Detector layout of the TA experiment. Filled 
black squares indicate the locations of the SDs, and trian¬ 
gles are the FD stations. The Central Laser Facility (CLF) 
is shown by the circle at the center of the array. 


5.1m^ spherical mirror, and the FD(BRM) and FD(LR) 
are newly designed telescopes with 6.8 m^ spherical mir¬ 
ror 0. Each FD station monitors the night sky from 
3° to 31° (MD) or 33° (BRM, LR) in the elevation, and 
112°(MD) or 108°(BRM, LR) in the azimuth. The read¬ 
out pixel of PMT camera is approximately 1.0° cone for 
all 3 stations. The PMT signal is recorded by the HiRes- 
1 sample and hold electronics at the MD station, and by 
a 10 MHz, 14-bit equivalent FADCs at the BRM and LR 
stations. A Central Laser Facility (CLF) is installed at 
the equidistant point from the 3 FD stations and regularly 
shoots a YAG laser into the sky for atmospheric monitor¬ 
ing. An Electron Light Source (ELS), a 40 MeV electron 
linac, is installed 100 m in front of the ED(BRM) to make 
an end-to-end calibration by shooting an electron beam 
into the sky. 

3 Energy Spectrum 

The E^ multiplied differential flux measured by the TA SD 
is plotted in Eigl2| where E is the energy of the primary 
cosmic rays. The spectrum is composed of 17,763 events 
above 10^^-^ eV collected in 6 years of SD operation from 
May 2008 to May 2014. The corresponding exposure is 
5400 km^ sr year above 10^^-^ eV. 

The analysis of the SD data is described in our previous 
publication |8| . Eor each SD event, the shower core loca¬ 
tion and the arrival direction were determined by fitting 
the shower lateral distribution and the SD hit timings with 
empirical formulae. Eor the spectrum measurement, the 
following cuts were applied to select events: 5 or more 
SDs, zenith angle 6 less than 45°, the core location more 
than 1200 m inside the array boarder, and the fitting 
and uncertainty estimates of energy and angle are within 
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Eigure 2: Energy spectrum obtained from the 6-year TA 
SD data. Solid black lines show the broken power law fit 
to the data. The red line is the extrapolation of the central 
segment to higher energies. 


certain limits O . The energy of the reconstructed event 
was determined by comparing the observed particle den¬ 
sity 800 m away from the shower axis, S(800), with the 
expected particle density estimated by the air shower sim¬ 
ulation. The CORSIKA and GEANT4 simulation tools 
and the de-thinning method were used for the primary 
proton to produce large numbers of air shower particles on 
the ground HU- The hadronic interaction model QGSJET- 
11-03 was used in the air shower generation. 

The obtained SD energy, Esd, is then compared with 
Efd, the energy measured by the ED for SD-ED simul¬ 
taneously measured (hybrid) events. The scatter plot of 
Esd/1-27 and Efd is shown in Eigl3]for energies above 
1018.5 Eigure, a simple rescaling of Esd 

restores a good equality between two measured energies. 

The aperture and the resolution of the accepted events 
were calculated by generating MC events according to the 
measured spectrum, and applying exactly the same anal¬ 
ysis and the selection criteria as the data. Distributions of 
simulated events for S(800), core location, energy, zenith 
and azimuth angles were checked for consistency with the 
distributions of data events. The acceptance is about 10% 
at 10^^-^ eV and rises to 100% above 10^^-^ eV. The en¬ 
ergy resolution above 10^^ eV is 20% and the directional 
resolution is 1.0° above 10^^-^ eV ll8ll^. 

As seen in Eigl2l the energy spectrum is well fitted by 
a broken power law (BPL) demonstrating the structure 
of “ankle” and “cutoff”. These structures were first ob¬ 
served by the HiRes experiment iTTOl in 2008. The PAO 
also saw a strong flux suppression around this energy in 
2008 ifTll . The fitted BPL parameters are listed in Ta- 
bl^ The energy E 1 / 2 , where the integral flux with cutoff 
becomes lower by a factor of 2 than the power-law extrap¬ 
olation, is 10^^-^^ eV for the data, and it agrees well with 
the prediction of the GZK cutoff (10^^-^^ eV) by Berezin¬ 
sky, Gazizov and Grigorieva 0. 
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Figure 3: Scatter plot of Esd/ 1-27 (abscissa) and Efd (or¬ 
dinate) for 551 hybrid events taken for May 2008 - May 
2013. The red line corresponds to the equality of the 
rescaled Esd and Efd- 


Table 1: Results of the BPL fit 


El 

1st break energy (ankle) 

j^Ql8.70+0.02 

E2 

2nd break energy (cutoff) 

19.74+0.04 

Ei/2 

GZK energy (1©) 

1019.73+0.04 ^Y 

n 

power index for E < Ei 

-3.30 + 0.03 

72 

power index for Ei < E < E 2 

-2.67 + 0.03 

73 

power index for E 2 < E 

-4.54 + 0.44 


Above the energy E 2 , 32 events are observed, whereas 
85.9 events are expected if the second break (cutoff) at 
E 2 does not exist. It signifies a 6.6 cr deviation from the 
continued spectrum without the cutoff. 

The observed spectrum is fitted with the expectation 
from the cosmic ray transport calculation and simulation 
ifT^ . The rate of cosmic ray production at redshift z is 
modeled as 

QiE, z) = aE-Pil+ E < E^,, z < z„ax (1) 

where p and m are the spectral index and the cosmological 
evolution parameter. The values of E^ax = 10^^ eV and 
Zmax = 2.0 were used such that these parameters do not af¬ 
fect the spectral shape in the fitting region. The generated 
protons are transported rectilinearly in the extra-galactic 
space, and energy losses caused by the interactions with 
the cosmic microwave background and infra-red photons 
were calculated. The calculation by the transport equation 
(TransportCR) and the MC simulation (CRPropa) were 
compared and they agreed well after some updates are ap¬ 
plied to the CRPropa parameters IT^ . 

The observed spectrum was fitted by optimizing 4 pa¬ 
rameters in the transport calculation: p, m, overall nor¬ 
malization (a) and the energy scale (AE). Eor the sim¬ 
ulated spectrum, events were generated, transported with 


Eigure 4: TA’s 5-year energy spectrum fitted with the sim¬ 
ulation of extra-galactic protons. Three cases of calcula¬ 
tion are shown: (Uniform) uniform cosmic ray source dis¬ 
tribution as described in the text, (ESS) source distribution 
according to the 2MASS Galaxy Redshift Catalog (XSCz) 
up to 250 Mpc, and uniform beyond 250 Mpc, (Uniform 
Zmin=0.01) uniform source distribution for 0.01 < z, but 
no source for z < 0.01. In plotting this Eigure, the best fit 
energy shift AE is not applied to the data, instead -AE is 
applied to the simulated spectrum. 


energy loss and the energy smearing same as the data reso¬ 
lution was applied before histogramming in the spectrum. 
Eor the data spectrum, only the statistical error was con¬ 
sidered in the fitting. The best fit to the data spectrum was 
obtained with p = 2.2, m = 6.4 and logAE" = - 0.05 

with x^NDE = 21.3/17. It is shown in EiglUas “Uni¬ 
form” B The best fit energy scale logAE" = - 0.05 is 

well within the systematic uncertainty of 21% for the TA 
hybrid events (H. 

FigH demonstrates that the TA 5-year SD spectrum 
can be fitted well with uniformly distributed, cosmolog- 
ically evolving cosmic ray proton sources without requir¬ 
ing contributions from other components. The allowed 
region of the p and m are shown in Eig|5l It should 
be noted that most particles above 10^^-^ eV are coming 
from z < 0.7, and constraints in Eig|5] apply only for such 
sources. 

4 Energy Spectrum by TALE 

The TA Low Energy extension, TALE, is a hybrid addi¬ 
tion to the TA installed next to the north ED station of TA, 
FD(MD) (see Fig©. It is composed of 10 additional ED 
telescopes refurbished from HiRes observing the higher 
elevation of 31° - 59° above the EoV of ED(MD). The 

^ The value of logA£ = -0.05 means the best fit was obtained when 
the measured energy is decreased by 10%. Adding the flux systematic 
error of ~3% in quadrature m changed the best fit to p = 2.2, m = 6.7 
and logA£ = - 0.03 with x^/NDF = 12.4/17. 
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Figure 7: TA and TALE energy spectrum with the BPL fit. 
TALE results in red open squares and in closed red circles 
are preliminary. HiRes data are from ifTOl . 


Figure 5: Permitted region of the source power index, p, 
and the evolution parameter, m, of uniformly distributed 
cosmic ray proton sources in z < 0.7. 


azimuthal coverage is ~100°. An infill array of 105 SDs 
is located in between the TALE FDs and the TA SD with 
400 m, 600 m and 1200 m spacings. The TALE SD works 
as an independent SD array as well as the hybrid array 
with TALE FDs. The TALE EDs were commissioned in 
September 2013. A total of 35 SDs were deployed and 
are being tested. The TALE and TA together will be mea¬ 
suring cosmic rays in a wide energy range of 10^^-^eV - 
1020.5 gY ^ single experiment. The use of Cherenkov 
analysis in TALE ED may decrease the energy threshold 
as low as 10^^-^ eV. 



Eigure 6: Detector layout of TALE. The TALE ED is in¬ 
stalled next to the ED(MD). The area between ED(MD) 
and the TA SD will be filled with in-fill arrays. 


In EiglTl the preliminary TALE energy spectrum is 


plotted together with TA SD, HiRes-1 and HiRes-2 mea¬ 
surements. The spectrum from 10^^-^ eV to 10^^-^ eV, des¬ 
ignated as TALE Ckov in Eig|7l is obtained by using the 
TALE ED as an Imaging Atmospheric Cherenkov Tele¬ 
scope (IACT) ifTbl . The Cherenkov events seen by TALE 
ED are with short duration (100 - 600 ns) and short track 
often confined in one ED camera. They are distinct from 
the fluorescence events, and are analyzed as a ED monoc¬ 
ular event. Its geometry is reconstructed by the profile 
constrained geometry fit (PCGE) developed for HiRes-1. 
The energy is obtained by using the Corsika/IACT simu¬ 
lation ifTTl . The spectrum from 10^^-^ eV to 10^^ eV are 
obtained by using the TALE ED as an air fiuorescence 
telescope. 

The whole spectrum in Eig|7]can be fitted by the BPL 
with two additional breaks at 10^^-^eV and 10^^-^ eV. Sim¬ 
ilar low energy features have been reported recently by 
KASCADE-Grande, Tibet ASy, IceTop, Tunka and other 
experiments (m. 

5 Composition 

Information on the primary cosmic ray composition is ob¬ 
tained from the measurement of Xmax in the atmosphere, 
which is expected to scale with InA and InE in average, 
where A is the mass number of the cosmic ray nuclide: 

< Xmax > — Cq, {InE — < InA >) -i- Cp (2) 

A most reliable way of measuring Xmax is the imag¬ 
ing of shower profile by the ED. Results so far obtained 
on Xmax are by the stereo ED measurement from HiRes 
d, the hybrid measurement from Auger (201, and the 
stereo measurement from TA EH. Each method corre¬ 
sponds to overlapping but specific phase space; the hybrid 
analysis requires the SD hits thus preferring more vertical 
events than the stereo, the stereo analysis requires higher 
energy threshold than the hybrid and so forth. Our goal is 
to establish a consistent set of results in the Xmax mea¬ 
surement among different analysis methods, and to make 
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Figure 8 : Xmax distribution for energies 10^^-^eV < E < 
10 i 8 . 6 eY Histograms are the results of MC simulations 
with QGSJET-II-03 hadronic interaction model: (Blue) 
iron and (Red) proton. MCs are normalized to the same 
total number of events with the data. Events are recon¬ 
structed as the hybrid of ED(MD) and SD. The data much 
more closely resembles proton simulation than iron simu¬ 
lation. 



Eigure 10: Elongation of average Xmax- Error bars are 
statistical only. MC expectation using several interaction 
models are shown for protons (blue) and irons (red). Both 
the data and the simulation include reconstruction and se¬ 
lection biases. 



Eigure 9: Same as Eigl 8 ]but for energies 10^^-^ eV < E < 
10 19.0 jYiQ data looks very protonic as in EigISj 


sure the differences between experiments are not due to 
the choice of different analysis methods. We report here 
the results of the hybrid analysis using the ED(MD) ll22ll . 

Eor this analysis, we selected events which were trig¬ 
gered and reconstructed independently by the monocular 
ED(MD) and by the SD. The ED and SD reconstructions 
are then combined into a hybrid analysis; the core location 
and impact timing information from the SD is combined 
with the timing information from the ED, thus improv¬ 
ing the geometry reconstruction of the monocular ED. 
The shower profile was then fitted with the Gaisser-Hillas 
function obtained from the pre-made CORSIKA shower 
library to determine the best fit Xmax- Pattern recogni¬ 
tion cuts were developed based on the experience of event 
eye scanning, and were used to remove events with poor 
fitting quality. 

We started with the good weather cut passed, fully re¬ 
constructed 1916 events for 10^^-^eV < E collected in 5 
years of TA runs between May 2008 and May 2013. After 
applying the geometrical cuts: bracketed Xmax, 0 < 56°, 
and the SD core location inside the array, 843 events re¬ 


mained. The final sample passing the pattern recognition 
cuts was 438 events. The Xmax resolution with pattern 
recognition cuts were studied by simulation. At 10^^-^ eV, 
the Xmax resolution was 35 g/cm^ without pattern recog¬ 
nition cuts, then improved to 23 g/cm^ with cuts. At 10^^-^ 
eV, the resolution was 15 g/cm^ 1 ^ . 

The distributions of Xmax for two energy regions E ~ 
10i8.5ev and E ~ 10^^-^eV are shown in Eig( 8 ]and Eig|9] 
together with the MC simulation results for protons and 
ironsEI. The MC events went through the same reconstruc¬ 
tion and selection procedures as the data. An elongation 
of Xmax is shown in Eig[T 0 | together with expectations 
for protons and irons using various hadronic interactions 
models: QGSJET-II-03, QGSJET-Olc and SYBILL 2.1. 

The measured average value of Xmax at 10 ^^ eV is 
751 + 16.3(sys) + 9.4(stat) g/cm^ and the elongation rate 
is 24.3 + 3.8(sys) + 6.5(stat) g/cm^. Assuming a purely 
protonic composition, taking into account all reconstruc¬ 
tion and acceptance biases and using the QGSJET-II-03 
model, we would expect the average Xmax at 10 ^^-^eV 
to be 763 g/cm^ and the elongation rate to be 29.7 g/cm^ 
per energy decade. We conclude that the measured com¬ 
position above 10 ^^-^eV is consistent with pure proton 
simulation: the same results as the stereo analysis using 
ED(BRM) and ED(LR). It is inconsistent with the pure 
iron simulation up to 10 ^^-^ eV, but above this energy the 
statistics are not enough to separate proton and iron. Con¬ 
sidering the uncertainty of the hadronic interaction mod¬ 
els, our results may not exclude the composition being 
light nucleus in this energy range. 

6 Anisotropy 

The PAO had shown in 2007 that arrival directions of 
UHECRs are correlated with AGNs in the VCV catalog 

^See (23 for other energy regions. 
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Figure 11: Red crosses show the time development of 
AGN correlation in the TA SD data. The abscissa is a cu¬ 
mulated number of events with energies 57 EeV < E and 
6 < 45.0°, the ordinate is the number of events correlated 
with the AGN within 3.1°. The AGNs with z < 0.018 
were taken in the VCV catalog as in the AGN correlation 
search by PAO. The blue line and regions represent the 
expectation and uncertainty from the PAO measurement 
121 . 



Figure 12: Sky plot of events with energies 57EeV < E 
and 6 < 55° in equatorial coordinates. The 5-year SD 
data is plotted. Galactic Plane (GP), Galactic Center (GC) 
and Super-Galactic Plane (SGP) are indicated. The bro¬ 
ken line shows the limit of TA acceptance at deck = -10°. 


in the southern hemisphere IH. The most significant 
correlation was obtained for events 60 EeV < E, associat¬ 
ing within 3.1° with the AGN less than 75 Mpc away (or 
z < 0.018). After the best fit condition was determined, an 
independent sample of events confirmed the chance prob¬ 
ability was 1.7 X 10“^. Eor the combined set of events 
with updated calibration, the largest correlation was ob¬ 
tained for events with 57 EeV < E. 

We repeated the correlation analysis with AGN in the 
northern hemisphere using the same selection criteria and 
the same catalog as the PAO measurement. Pig [TT] shows 
the chronological progress of the number of correlated 
events collected in 5 years of TA SD runs. As of May 
2013, we obtained 42 events above 57 EeV and 17 events 
were correlated with the AGNs within 75 Mpc whereas 
we expect ~10 events are accidentally associated with 472 
AGNs in the catalog if the events are randomly distributed 



Pigure 13: Skyplot of events in the 6th year only 



Pigure 14: Significance of fiux enhancement. Oversam¬ 
pling with 20° radius circle is made for events in Pig US 
The random background was calculated by throwing 100k 
sets of 72 isotropic MC events. The maximum signifi¬ 
cance is reached at R.A.=146.7° and decl.= 43.2°. The 
Li-Ma significance of this spot is 5.1 cr. 


over the sky. The probability of 17 or more events corre¬ 
lated by chance was 1.4% ll24l . This analysis used the 
same sample as the SD spectrum analysis with tight cuts: 
6 < 45.0° and the shower core to be more than 1200 m 
inside the array border. 

In order to quantify the significance of a possible fiux 
enhancement seen in the highest energy sample of 57 EeV 
< E off the Super Galactic Plane |(25l, we loosened the 
event selection to accept all events with 0 < 55° and with 
the core location inside the array border. It increased the 
sample size above 57 EeV from 42 events to 72 events in 
the 5-year sample for the “hotspot” analysis (Pig US. The 
looser cuts deteriorated the angular resolution from 1.0° 
to 1.7°; the energy resolution became from 15% to 20% 
according to the simulation. 

We made an oversampling of the 72 events with a cir¬ 
cle of radius 20.0° in the sky plot. An average of expected 
number of events from the isotropic arrival direction was 
calculated for each point of the sky by producing 100k 
sets of 72 MC events uniformly thrown over the sky (ex¬ 
posure oc sin(^) cos(^)). The maximum deviation from the 
isotropic expectation was then searched over the sky with 
0.1° step. 

The obtained significance map is shown in Pig [141 The 
largest deviation of 5.1 cr occurred at R.A.=146.7° and 
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decL= 43.2°, where 19 data events were found within the 
circle of r = 20°, and 4.49 events were expected from the 
isotropic arrival distribution ll^ . 

We estimated a chance probability of this “hotspot” to 
appear from the isotropic arrival distribution by simula¬ 
tion: IM sets of 72 events with isotropic distribution were 
generated and analyzed in the same way as the data by 
searching for the sky direction to give the maximum de¬ 
viation. In the search for the maximum deviation, we not 
only searched with r=20° but also included the searches 
with r = 15°, 25°, 30°, 35°. This is to compensate a poten¬ 
tial “eye-ball scan” we might have made unintentionally 
in selecting the radius of 20°. 

Out of IM sets simulated, we found 37 sets whose max¬ 
imum significance exceeded 5.1 cr. The chance probabil¬ 
ity of the hotspot appearing from the isotropic distribution 
somewhere in the TA’s acceptance (direction of hotspot) 
and with the radius of 15°- 35° (size of hotspot) was de¬ 
termined to be 3.7xlO""^ or 3 Act. 

We repeated the same analysis for the data collected in 
the 6th year. May 2013 - May 2014. We obtained 15 ad¬ 
ditional events above 57 EeV by the same loose cuts (see 
FigfT3l) and 4 events are found in the hot spot determined 
by the 5-year data sample. For all the data collected in 6 
years, the significance of the hotspot increased to 5.55(T 
and the chance probability for getting more than 5.55(T 
from the isotropy advanced to 4.0 cr (from the 3.4 cr of 
the 5-year data). 


7 Prospects 

TA plans to extend the acceptance at the highest energy 
region by adding 500 additional SDs and 1 or 2 more FD 
stations. The plan is called TAx4 and the detector layout 
is shown in Fig [151 The SDs of the same design will be 
deployed with a larger spacing of 2.1km extending east of 
the existing TA SD array. The new array will be 96% effi¬ 
cient at 57 EeV for the trigger and reconstruction, and the 
SD acceptance will be quadrupled to -3000 km^ above 
this energy. The efficiency at 10^^ eV is expected to be 
26%. The new FD stations will be composed of refur¬ 
bished HiRes telescopes with FADC readout already used 
for the TAFE FD telescopes. 

The construction of TAx4 will take 2 years, and the 
corresponding budgets have been applied in Japan, USA, 
Korea and Russia. After 3 years of operation in 2018- 
2020, TAx4 will accumulate 21-year equivalent of TA SD 
data and 18-year equivalent of TA hybrid data, if the con¬ 
struction is approved in 2015. With -250 events collected 
above 57 EeV by 2020, the hotspot may be established 
at 7 (T or higher and properties of the hotspot such as the 
energy spectrum, the composition and finer sub-structure 
will be studied. There may be other weaker hotspots ap¬ 
pearing in the northern sky. The Xmax information will 
be available in the cutoff region, and the cutoff shape will 
be studied in detail with respect to the expectation of the 
GZK cutoff Il27l. 



Figure 15: Planned detector layout of TAx4. Additional 
SDs and FDs will be installed on the east side of the ex¬ 
isting TA. The SD area will be quadrupled for the highest 
energy events. 


8 Summary 

In this conference, we presented the energy spectrum and 
the anisotropy measured by TA using 6 years of SD data. 
We also presented Xmax results analyzed in the hybrid 
mode. 

The observed spectrum above 10^^-^eV shows a dip at 
IQis.VeV and a cutoff at 10^^-^eV. These features are well 
reproduced by the model calculation of the GZK process: 
energy loss of extra-galactic protons by the interaction 
with the CMB and IR background. The simulation sug¬ 
gested a source spectrum of and source evolution of 
(1-i-z)^-^ for z < 0.7. 

A comparison of Xmax distribution with model simula¬ 
tions (QGSJET-II-03), we showed the primary composi¬ 
tion is consistent with 100% proton and inconsistent with 
100% iron for energies 10^^-^ eV < 10^^-^ eV. Due to the 
ambiguities arising from hadronic interaction models, we 
can not refute an inclusion of light nuclei in the same en¬ 
ergy range. The energy region above 10^^-^ requires more 
statistics for a better suggestion of the composition. 

Above the cutoff energy, 57EeV < E, we observed a 
hotspot near the direction of Ursa-Major with a size of 
r=20°. The probability of this hotspot being formed by 
chance from the isotropic distribution is 4 cr. The center 
of the hotspot is 19° off from the Super-Galactic Plane, 
and there is no immediate candidate of the UHECR source 
nearby. 

Further studies of UHECRs at TA will continue with an 
extension project TAx4, which will quadruple the SD ac¬ 
ceptance and double the hybrid acceptance at the highest 
energy by 2017. 
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